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We perform simulations of (3-colour) QCD with 2 quark flavours at a finite chemical 
potential /ij for isospin(7a), and of 2-colour QCD at a finite chemical potential fi for quark 
^"*^ I number. At zero temperature, QCD at finite fii has a mean-field phase transition at fii — 

to a superfluid state with a charged pion condensate which spontaneously breaks ^3. We 
study the finite temperature transition as a function of fii. For /17 < m-x, where this is 
closely related to the transition at finite /i, this appears to be a crossover independent of 
quark mass, with no sign of the proposed critical endpoint. For fii > this becomes 
a true phase transition where the pion condensate evaporates. For fii just above m w the 
transition seems to be second order, while for larger fii it appears to become first order. At 
zero temperature, 2-colour QCD also possesses a superfluid state with a diquark condensate. 
We study its spectrum of Goldstone and pseudo-Goldstone bosons associated with chiral and 
quark-number symmetry breaking. 
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§1. Introduction 

We are interested in QCD at finite baryon-number, isospin and strangeness den- 
On. sities at zero and at finite temperature. Finite baryon-number is of most interest, 
but most difficult to simulate. QCD at finite chemical potential fx for quark number 
has a complex fermion determinant which prevents the use of standard simulation 
methods. 

We are studying theories which possess some of the properties of QCD at finite 
fi, but have positive fermion determinants. This allows us to use standard (hybrid 
molecular-dynamics) simulations. 

2-colour QCD at finite \x has a positive fermion pfaffian and is amenable to sim- 
ulations.^^ It has a sensible meson spectrum, but unphysical 'baryon' spectrum. 
At fi = it undergoes a phase transition to a superfluid state with a diquark 

condensate and Goldstone bosons P^^^ This is the analogue of the colour super- 
conducting state suggested for QCD at finite fi. We have also studied this theory 
at non-zero temperatureP Here we summarize our calculation of the spectrum of 
(pseudo)-Goldstone bosons which characterize the pattern of symmetry breaking for 
this theory™ 

QCD at a finite chemical potential [i\ for isospin (I3), has a positive fermion 
determinant which allows simulations .™ It undergoes a phase transition at [ij = 
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to a superfluid state with a charged pion condensate which breaks 1% and parity 
spontaneously, with associated charged Goldstone pionsP'^ 

Let us now compare the proposed phase diagrams for QCD at finite /JJHO and 
QCD at finite [ij (figure Both show crossovers (dashed lines) at small chemical 
potential. For QCD at finite /x ^ the 2 solid lines which intersect the T = axis 
and end in open circles (critical points) are expected to be first order. The order of 
the transition between the colour-superconducting and quark-gluon plasma phases 
is unknown. For QCD at finite fij, the solid lines to the left of the open circle 
(tricritical point) are second order, while that to the right is first order. 



N f = 2 QCD at finite baryon density N f =2 QCD at finite isospin density 




m N /3 fj, m„ /Uj 

Fig. 1. a) Proposed phase diagram for QCD at finite quark-number chemical potential p and tem- 
perature T. b)Proposed phase diagram for QCD at finite isospin chemical potential /ij and 
temperature T. 



We are studying 2-flavour QCD at finite m and temperature in the neighbour- 
hood of the finite temperature transition.™*^ At low [ij this transition from 
hadronic matter to a quark-gluon plasma is related to the corresponding transition 
at finite [i. Our predictions for the small /i behaviour of the finite temperature tran- 
sition can be compared directly with those obtained by de Forcrand and Philipsen 
by continuing from imaginary //P^J These results will also be compared with results 
obtained using the series methods of the Bielefeld-Swansea collaborationPJ'EJ They 
could also be checked using the reweighting methods of Fodor and KatzP^*^*^3 
We are also simulating this theory at large \ii where the finite temperature transi- 
tion occurs at the point where the pion condensate evaporates. It is therefore a true 
phase transition. In particular we wish to determine where this transition changes 
from second order to first P 

Section 2 gives our results for lattice QCD at finite fii and zero temperature. 
In Section 3 we present preliminary results for QCD at finite fii and temperature. 
The (pseudo)-Goldstone spectrum for 2-colour lattice QCD at finite fi is given in 
section 4. In section 5 we indicate how our new xQCD action might help overcome 
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some of the problems at finite h/m- Section 6 gives our conclusions and suggests 
future directions for investigation. 

§2. Lattice QCD at finite /xj 



The staggered quark action for lattice QCD at finite \ii is 



E 

sites 



(2-1) 



where L\^t^lii) is the standard staggered quark ^with the links in the +t direction 
multiplied by exp^r^Lij) and those in the —t direction multiplied by exp(— |t3//j). 
The A term is an explicit 1$ symmetry breaking term required to see spontaneous 
symmetry breaking on a finite lattice. The determinant 

det[#(-73/i/) + m + iXer 2 ] = det[,4M + A 2 ], (2-2) 

where 

A = fl(-»i) + m, (2-3) 

is positive allowing us to use standard hybrid molecular dynamics simulations. 

This lattice action has a global £7(2) x U(2) flavour symmetry when m = = 
A = 0. This breaks spontaneously to U(2), with 4 Goldstone pions. When lii = but 
this is reduced to U(2)y, and the would-be-Goldstone pions gain a mass by 
PCAC. If m = 0, but in / 0, the symmetry is reduced to U(l) x U(l) x U(l) x U(l). 
Finally, if m ^ and \ii 7^ 0, the symmetry is further reduced to U(l)y x U(l)v- 

As [ii is increased from zero, the effective mass of the tt + is reduced to 

m^m) = m.* - m, (2-4) 

which vanishes at [ii = m n . For lij > m n , the U(l)y symmetry associated with T3 
is broken spontaneously by a charged pion condensate 

i(x^2X) O i(^l5T2^)- (2-5) 

The charged pion excitation created by 

ixerix (2-6) 

is the associated Goldstone boson. When A is non-zero this gains a mass. Note that 
if m = 0, there will actually be 2 Goldstone bosons. 

Simulations of lattice QCD at finite lli were performed with m = 0.025 and A = 
0.0025, 0.005 on an 8 4 lattice.^ Figure |21 shows the diquark and chiral condensates as 
functions of lij. For m n and \xj small one can use chiral perturbation theory to predict 
the behaviour of this theory At tree-level this shows the phase transition at 
lii = m n . For lii > m n , the condensate rotates from the chiral direction towards the 
pion condensate direction so that 1% breaks spontaneously. Indications are that at 
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Fig. 2. a) Charged pion condensate as a function of ^3 chemical potential fii. The curves are a fit 
to mean- field scaling, b) Chiral condensate as a function of chemical potential/!/ . The curves 
are predictions from the fits to the pion condensate. 



higher orders, the transition remains at fii = m n , but the rotation of the condensate 
is accompanied by an increase of its magnitude. However, the transition remains 
mean-field. We model this behaviour by fitting to a scaling form based on the tree- 
level predictions of a chiral Lagrangian of the linear sigma model form which allows 
the magnitude of the condensate to vary, while retaining the other predictions of 
chiral perturbation theory, with only one additional fitting parameter. The fits are 
shown on the figures. These indicate that there is a phase transition with mean-field 
exponents to a state where I3 and parity are broken spontaneously. These fits also 
predict how the ^3 density increases from zero above m n , and are in accord with the 
data (not shown) provided /ij is not too large. 

§3. Lattice QCD at finite /ij and T 

We first consider the finite temperature transition for /u/ < \x c . As fij is increased 
from zero, the (3 C of the crossover from hadronic matter to a quark gluon plasma 
decreases. Note, the results in this section should be considered preliminary. 

We have performed simulations on an 8 3 x 4 lattice at m = 0.05 and A = 0, at 
values of fij up to 0.55, which is very close to /j,J^^ The position of the crossover 
is determined from the maxima of the susceptibilities for the plaquette, the chiral 
condensate, the Wilson line (Polyakov Loop) and the isospin density. As we shall see, 
there is consistency between these estimates. (The close proximity of the transitions 
in these various observables has been noted by many authors. Recently attempts 
to understand the reason for this have been made for zero /i///jP^ Combining this 
with earlier observations for finite //j^ it seems probable that this argument works 
at finite n (and hence fj,j). A detailed analysis in terms of effective field theories, 
which also indicates applicability to finite has been given iiP^^ ,) "w e use 
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Ferrenberg-Swendsen reweighting**' to determine these values with some precision. 
For the chiral condensate and isospin density which are determined from stochastic 
estimators, we use 5 noise vectors to calculate the susceptibilities, discarding the 
diagonal terms to obtain an unbiased estimate. 

In figure 01 we show the Thermal Wilson Line (Polyakov Loop) as a function of 
(3 = 6/g 2 for a selection of \ij values in this range. First is is clear that for each 
Hi there is a rapid crossover from hadronic matter to a quark-gluon plasma over an 
interval A/3 ~ 0.05. What we also notice is that the (3 of the crossover decreases 
monotonically with increasing 



SU(3) N,=2 m=0.05 A=0 8 3 x4 lattice 



SU(3) N,=2 m=0.05 X=0 8 3 x4 lattice 
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Fig. 3. Wilson Line as a function of /3 = 6/g 2 
for /x/ values < 0.55. 



Fig. 4. Isospin density as a function of P = 
6/g 2 for /x/ values < 0.55. 



The corresponding graph for the ^3 density is given in figure |3J Here, not only 
does the (3 value for the cross over decrease with increasing m, the value of this 
density increases with increasing /zj. This is expected, since when m is increased, 
it raises the level of the fermi sea, and hence the excess of ^3 = 1/2 quarks over 
I3 = —1/2 quarks. 

In figure we show the Wilson Line susceptibility and in figure we show the 
susceptibility for the isospin density. The chiral and plaquette susceptibilities are 
similar to that for the Wilson Line. The susceptibility for an operator O is defined 

as 

Xo = V{0 2 - (O) 2 ), (3-1) 

where V is the space-time volume of the lattice. The peak of the susceptibility 
serves to define the position of the crossover. Clearly the measurements in these 
plots are inadequate to give precise positions of the peaks. However we can make a 
Ferrenberg-Swendsen reweighting of our 'data' at (3 values close to the peaks, which 
yields a more precise estimate of their positions. 

We plot our estimates of (3 C versus /x| from these 4 susceptibilities in figure UJ 
since we expect the leading dependence to be quadratic in \xj. The straight line in 
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Fig. 5. Wilson Line susceptibility as a func- Fig. 6. Isospin susceptibility as a function of 



tion of (3 = 6/g for fii values < 0.55. 



j3 = 6/g for fii values < 0.55. 
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Fig. 7. /3 C determined from the maxima of the susceptibilities, as a function of fij . 
this graph 

& = 5.322 -0.143/x?, (3-2) 
is not a proper fit and is only meant to be a rough guide. If we assume that at 
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fj,j = 0, T c « 173 MeV, then 2-loop running of the coupling would predict that for 
the highest \xi considered, T c (/xj = 0.55) ~ 164 MeV. Hence the relatively weak 
dependence of (5 C on m does translate into a relatively small change in T c over a 
fairly large range of hi. 

The Swansea-Bielefeld collaboration have determined that, for quark-number 
chemical potential /j, small enough, the phase of the fermion determinant is well 
enough behaved that simulating with the magnitude of the determinant and includ- 
ing the phase factor in the measurements should work. This means that, to the 
extent that the position of the crossover is a well-defined quantity (i.e. is the same 
for all observables) , the dependence of T c on fi and fii should be identical for small 
fi, in particular, 

= fl^j = lii) (3-3) 

so we can predict the small ji behaviour of (5 C and hence T c . Our results in equa- 
tion 13-21 using equation 13-31 appear to be in agreement with those of de Forcrand 
and PhilipserP2» 

This also suggests that we should look for a critical point where the crossover 
becomes a first order transition, since such a critical endpoint is predicted for QCD 
at a finite chemical potential fi for quark number. For the m = 0.05 simulations 
reported here, we see no evidence for this for < 0.55 (i.e. for fij < /i c ). However, 
Hi = 0.55 is equivalent to \i = 0.275 rs 180 MeV {^B = 535 MeV), which is less than 
H measured for the critical point by Fodor and KatzP^^^ 1 ^ (Note, however, that 
the Fodor-Katz simulations are for 2+1-flavours, not 2.) We are therefore performing 
simulations with m = 0.1 and m = 0.2 where \x c is beyond the Fodor-Katz value. 
However, we still see no signals (so far) for first-order behaviour. Thus we conclude 
that if there is a transition to be identified with this critical point, its m increases 
with increasing mass, which is not unexpected. 

Now let us turn to our studies of the finite temperature transition for larger fj,js. 
For hi > He, the pion condensate evaporates at the finite temperature transition. 
Thus here the finite temperature transition is a phase transition. At low tempera- 
tures, by approaching this phase boundary in the fil direction we have determined 
that the transition is second order. We now have preliminary evidence that for 
m = 0.05, A = 0.005 and Nt = 4, the fii = 0.8 transition is first order, based on 
simulations on 8 3 x 4 and 16 3 x 4 lattices, and are in the process of determining the 
position of the tricritical point. 

We note that this transition is at (3 C ~ 5.2675 i.e. at T c ~ 162 MeV. This \jli 
value is equivalent to fi = 0.4 ~ 259 MeV. Hence the tricritical point must be close 
to the Fodor-Katz critical endpoint where T c ~ 160 MeV, [i ~ 242 MeV. We can 
therefore speculate that these 2 transitions are somehow related. Figure |HK shows 
the pion condensate function of (3 on a 16 3 x 4 lattice at fi = 0.8. Figure 
shows the time evolution of this condensate from hot and cold starts at (3 = 5.2675, 
showing evidence for coexisting states. Our "time" increment for updating dt = 0.02 
was larger than desirable so we are repeating our simulations with dt = 0.01 close to 
the transition. 
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a TIME 

Fig. 8. a) Pion condensate as a function of j3 — 6/g 2 for m = 0.05 and /ij = 0.8. b) Pion condensates 
as functions of molecular-dynamics time at f3 = 5.2675 from hot and cold starts. 



§4. The Spectrum of 2-colour QCD 

Lattice 2-colour QCD at finite chemical potential [/, for quark number has the 
staggered fermion action 

S f=Yl WW + m \x + \KX T T2X + XT 2 x T }\ (4-1) 

sites 

This theory has a positive pfaffian which allows us to simulate it using the hybrid 
molecular-dynamics method. This lattice action has a U (2) flavour symmetry when 
m = fi = A = 0. When this breaks spontaneously U{2) — ► U(l) giving rise to 
3 Goldstone bosons. We have studied the behaviour of these bosons when these 
parameters are no longer zero. In particular we study the \i dependence when m is 
fixed at a value small enough to be within reach of chiral perturbation theory, and 

a « m jam 

At A = this theory has a phase transition at fj, = to a superfluid phase 

with a diquark condensate which spontaneously breaks quark- number. There is 
1 true Goldstone boson. As for QCD at finite ///, the condensates and quark- 
number density are well described by a tree-level Lagrangian of the linear sigma 
model class. Figure El shows the diquark condensate, which is an order parameter 
for this transition as a function of \i and A, as measured in simulations on a 12 3 x 24 
lattice at = A/g 2 = 1.5 and m = 0.025. 

We measure the \i and A dependence of the 3 pseudo-Goldstone bosons and 
compare them with the predictions of the linear sigma model effective Lagrangian. 
{m-K is small enough at m = 0.025 for us to expect chiral perturbation theory to 
be valid.) Figure ITOl shows the 3 (pseudo)-Goldstone boson masses as functions of 
[i. Although there is good qualitative agreement, quantitative agreement is lacking 
once // is appreciably greater than m n /2. 
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SU(2) N,=4 /S=1.5 m=0.025 12 3 x24 lattice 
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Fig. 9. Diquark condensate as a function of \x. 
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§5. xQCD at finite and T. 

We have formulated QCD with additional irrelevant chiral 4-fermion interactions 
which make the Dirac operator non-sing ular in the chiral limit (xQCD)PJ The 
staggered fermion part of this action is: 



S f 



sites L i ) s 



where s are the sites on the dual lattice, and the sum over i is the sum over the 16 
sites on the dual lattice which are 1 unit away from the given site on the original 
lattice. It is just the properties which allow this action to be used in the chiral limit, 
that promise to make it better behaved than the standard staggered action at finite 
chemical potential. 

Now consider what happens when we add a chemical potential fx in the usual 
way. Since the quark propagator on a single configuration no longer behaves as 
exp[— (m 7r /2)\x — y\], but is controlled by ((a) = (xx)/l) + m rather than m, the 
phase of the determinant can remain well behaved for [i > m n /2, and even in the 
chiral limit. 

Thus, provided we keep 7 from becoming too large, the Swansea-Bielefeld ex- 
pansion should be well behaved, even for m = 0. The hope is that the phase of the 
determinant will remain under control up to the critical end point. Since keeping 7 
"small" makes the flavour symmetry breaking bad, this should be considered as a 
model. It remains to see how well it represents QCD at finite //. 

In addition, we should be able to use this action to simulate with the magnitude 
of the determinant (equivalent to simulating at finite hi) at m = and to search 
for the tricritical point. Here we see that the pion consisting of a quark from the 
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Dirac operator and an antiquark from the conjugate Dirac operator will be relatively 
heavy, which prevents it from causing problems at small /x. 

Preliminary studies of this theory by Barbour, Kogut and MorrisorP^ suggest 
that xQCD is indeed better behaved at finite /x than standard lattice QCD, and are 
consistent with the scenario presented above. 

§6. Conclusions 

QCD at finite isospin chemical potential and 2-colour QCD at finite quark- 
number chemical potential have some of the properties of QCD at finite quark/baryon- 
number chemical potential. Unlike QCD at finite /x, they have positive fermion 
determinants, so they are amenable to standard simulation methods. 

QCD at finite fii has a phase transition with mean-field critical exponents to 
a superfluid phase in which isospin (J3) and parity are broken by a charged pion 
condensate at ti/ = m n , as predicted by chiral perturbation theory. The dependence 
of the position of the finite temperature transition (T c ) on ti and on fij is the same 
for ii sufficiently small; so simulations at finite /x/ are a convenient way of estimating 
the [i dependence. We will check the agreement between our measurements and the 
Swansea-Bielefeld predictions. It appears that our results are consistent with those 
obtained by de Forcrand and Philipsen from simulations at imaginary xi. 

Our 2-flavour simulations do not show any sign of the critical endpoint expected 
for QCD at finite /x and hence /i/. Perhaps the tricritical point on the boundary of 
the pion condensed phase, where the transition changes from second order to first 
order, is a remnant of this endpoint. Unfortunately this is in the region where we 
cannot relate the finite /x/ and finite /i behaviour of lattice QCD. 

For 2+1 flavours Fodor and Katz were able to determine the position of the 
critical endpoint at finite ^GSEDEfl De Forcrand and Philipsen have made a deter- 
mination of the position of this critical endpoint for 3 flavours by continuation from 
imaginary /xP^ (Similar determinations have been made for 4 flavours by D'Elia 
and Lombardo. 29 ') The Bielefeld-Swansea collaboration have also determined the 
position of the critical endpoint for 3 flavours using their series expansions around 
/x = Op" In addition to the fact that 2+1 flavours describes the real world, 2+1 or 3 
flavours have the advantage that for small enough quark masses, the critical endpoint 
is at the origin. This means that one can arrange for the critical endpoint to be as 
close to [i = as one desires by tuning the quark masses, and follow its trajectory as 
these masses are increased. In particular, in a 3-flavour version of QCD at finite /xj, 
we should be able to arrange for the critical endpoint to occur for /x/ < m w , where 
it should be related to the critical endpoint at finite /x. We are currently extending 
our QCD at finite /x/ simulations to 3 flavours, ignoring the fact that the theory 
we are simulating is not expected to have a sensible continuum limit, and hope to 
determine the position of this endpoint with good precision. 

There remain some questions to be answered concerning QCD at finite /x/. These 
include: Does pion condensation at finite isospin density occur when the system is 
also at finite baryon-number density? If so, how large can the baryon-number density 
become before pion condensation no longer occurs? How well does the pseudo- 
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Goldstone spectrum of lattice QCD at finite [ii compare with the predictions of 
chiral perturbation theory? 

The pseudo-Goldstone spectrum of 2-colour lattice QCD at finite [i shows semi- 
quantitative agreement with the predictions of effective chiral Lagrangians. 

The xQCD action shows promise for simulations at finite [i at T . It should allow 
determination of the critical endpoint for small (or even zero) quark masses. 
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